Slow dynamics in a driven two-lane particle system 
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We study a two-lane model of two-species of particles that perform biased diffusion. Extensive 
numerical simulations show that when bias q is strong enough oppositely drifting particles form 
some clusters that block each other. Coarsening of such clusters is very slow and their size increases 
logarithmically in time. For smaller q particles collapse essentially on a single cluster whose size 
seems to diverge at a certain value of q = q c . Simulations show that despite slow coarsening, the 
model has rather large power-law cooling-rate effects, ft makes its dynamics different from glassy 
systems, but similar to some three-dimensional Ising-type models (gonihedric models). 
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Despite its apparent simplicity, statistical mechanics of 
low-dimensional driven diffusive particle systems is still 
attracting a considerable interest This is due to nu- 
merous applications ranging from vehicular or pedestrian 
transport Q , to gel electrophoresis [U , to molecular mo- 
tors [J]. Moreover, such systems exhibit the wealth of 
highly nontrivial and often surprising features. Indeed, 
in the steady state of these systems we observe novel uni- 
versality classes [|| , off-critical long-range correlations Q , 
or internal energy that is a decreasing function of tem- 
perature Q- It shows that experience we have gained 
studying equilibrium systems is not very helpful in the 
realm of nonequilibrium phenomena. 

There is a growing evidence that dynamics of these 
nonequilibrium systems is also very interesting. For ex- 
ample, asymptotic cluster growth in a certain two-species 
model is governed by the exponent that is twice as large 
as could be (naively) predicted using Lifshitz-Slyosov the- 
ory 0] . Nucleation of clusters in this model is also sur- 
prising: residence distribution develops a peak that sig- 
nals formation of a macroscopic cluster, but that peak 
disappears, however, when the system size becomes ex- 
tremely large @. 

For some other models Evans et al. predicted that 
power-law coarsening, that typically accompanies sym- 
metry breaking transitions, should become logarithmmi- 
cally slow. However, due to numerical difficulties such a 
behaviour was observed only in an effective toy model. 
Such a slow coarsening is one of the important ingredients 
of glassy dynamics. Absence of quenched disorder means 
that kinetic barriers which slow down the evolution are 
generated by the dynamics of the model. A similar sce- 
nario takes place in the so-called ordinary glasses pd| . 
Another important ingredient of glassy dynamics are very 
small cooling-rate effects, i.e., very slow (presumably log- 
arithmically slow) growth of the characteristic length 
scale as a function of the inverse of the cooling rate. 
Some three-dimensional Ising-type models exhibit such 
a behaviour [13]. Although generation of barriers (even 
diverging ones) might lead to slow-coarsenin g it does not 
necessarily imply small cooling-rate effects [13| . In ad- 
dition, some arguments suggest that in two-dimensional 



Ising-type models thermal fluctuations should suppress 
generation of energy barriers restoring fast, power-law 
coarsening [l4| . The role of dimensionality in the glassy 
dynamics is important also in off-lattice systems. Indeed, 
recent Molecular Dynamics simulations of one-, two- 
, three-, and four-dimensional Lennard- Jones systems 
show that the tendency to form glasses increases with di- 
mensionality and in one-dimensional systems glassy tran- 
sition was absent (l5j . Let us notice, however, that sim- 
ulations of off-lattice systems are computationally very 
demanding and it is often mere presence or absence of 
hysteresis that is used to identify the glassy transition. 

The main objective of the present paper is to check 
whether absence of glassy transition in one-dimensional 
systems might be a property also of lattice models. 
Slowly coarsening driven diffusive particle systems are 
good candidates for such an analysis. To examine 
whether other ingredients of glassy dynamics might ap- 
pear in such low-dimensional systems it would be de- 
sirable to analyse appropriately defined cooling rate ef- 
fects. Since analytically tractable driven diffusive mod- 
els usually can be solved only with respect to steady 
state properties [16j |. further understanding of dynam- 
ics of these systems most likely have to rely on numerical 
approaches. In systems with slow dynamics development 
of efficient computational methods is particularly impor- 
tant. In the present paper, we examine a model where 
two species of particles perform a biased diffusion on a 
one-dimensional (two-lane) lattice. Closely related mod- 
els have already been analysed in the context of formation 
of spatial structures [l7j . nucleation kinetics, and cluster 
growth [1, Q and our work contributes to better under- 
standing of this class of models. Our simulations with an 
efficient algorithm clearly indicate that the regime with 
spatial inhomogeneities is characterised by logarithmi- 
cally slow coarsening. However, rather large, power-law 
cooling-rate effects show that phase transition separating 
homogeneous and inhomogeneous phases of the model 
differs from realistic glassy transition. We also discuss 
similarities of the dynamics of our model and that of 
some Ising-type models. 

In our model there are N particles of two kinds, 
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FIG. 1: Probabilities of elementary transitions in our model. 
Forp2 > pi positive/negative particles are preferentially mov- 
ing to the right /left and some gridlocks might form as the 
one to the right. However, collective movements of particles 
might create some holes and eventually eliminate such grid- 
locks. Since such movements are against the drive and thus 
are highly unlikely, the lifetime of such gridlocks rapidly in- 
creases with their size. 

" positive" and " negative" , that perform biased diffusion 
on a two-lane lattice of size 2xL (all results reported 
in this paper are obtained for density = 0.1). Each 
site of the lattice can be occupied by at most one 
particle, and we assume that the number of "positive" 
and "negative" particles are equal. In addition, particles 
are exposed to "electric" field that introduces some bias 
on their diffusion in the horizontal direction. Detailed 
dynamics is specified below: 

(i) Select randomly a particle. 

(ii) Select the target site (one of the three nearest 
neighbours of the selected particle) with probability 
depending on the kind of the selected particle: For 
" positive" particles left and right neighbours are selected 
with probabilities p\ andp2, respectively; For "negative" 
particles these probabilities are exchanged (see Fig[TJ. 
The probability of inter-lane jumps p^ is for both kinds 
of particles the same and all results described in the 
present paper are obtained for p^ = 0.5. An important 
parameter of the model, that is related with the strength 
of the "electric" field, is the bias probability q defined 
using the difference of probabilities of right and left 
jumps as q — |(p 2 ~Pi), and of course, p\ +P2 +Pa = 1- 
Let us also notice that dynamics preserves CP symmetry 
(i.e., probability of jumping to the right for a "positive" 
particle is the same as that of jumping to the left for a 
"negative" particle). 

(iii) Provided that the target site is empty, the selected 
particle jumps to that site. Otherwise, the attempt to 
move the selected particle is rejected. 

Closely related models have already been analysed by 
Schmittmann et al. [H, HI, Il7| . In their formulation of 
the dynamical rules, however, it is a lattice bond that is 
selected randomly. As a result, their algorithm is efficient 
but only for dense systems (when density is small there 
is a large probability of selecting a bond with no particles 
at its ends and no moves are possible) . In our case we can 
ensure that each time we select a particle. However, when 
large clusters of particles are formed, a lot of particles 
get blocked. In such a case it is useful to keep the list 
of active particles (i.e., those that can move), and select 
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FIG. 2: The time evolution of a random initial configuration 
of TV = 5 ■ 10 s particles on a ladder of length L — 25 ■ 10 3 and 
for q = 0.045. After around t = 5 ■ 10 5 MC steps essentially 
all particles form a single cluster. The configuration of every 
fifth particle is recorded every 7500 MC steps. 

particles only from such a list [l^|. We have to notice, 
that since a unit of time in our model is defined as a 
single, on average, update of each particle, selecting only 
from the list of active particles the simulation time t is 
increased by 1/Na, where Na is the (current) number 
of active particles. Our algorithm was actually a hybrid 
version: when the number of active particles is large, 
we select a particle out of all particles, and only when it 
drops below the threshold we switch to the version where 
we select out of the list of active particles (as a threshold 
we used the value Na = N/5). Let us also notice that in 
the late-stage evolution there is usually a small number 
of active particles and the described algorithm is much 
more efficient than that where selection is made out of 
all particles. 

It is already known [TtJ that for models of this kind for 
sufficiently small but positive bias q the model remains 
in the homogeneous phase with a slow drift of particles 
("positive" to the right and "negative" to the left). For 
larger q it was reported [l~7| that after some time a large 
cluster of particles is formed and the drift stops. Forma- 
tion of such a cluster is illustrated in Fig. [5]that shows the 
time dependence of the initially random configuration of 
5 • 10 3 particles evolving at q = 0.045. Indeed, one can 
notice that around t = 2.5 ■ 10 5 a cluster is formed that 
grow and after t ~ 5 • 10 5 almost all particles collapse on 
that cluster. 

Main results reported in the present paper are inspired 
by the observation that for larger bias q instead of form- 
ing a single cluster, the model gets trapped in the slowly- 
coarsening state with the wide spectrum of clusters of 
various sizes. An example of such an evolution is shown 
in Fig. [3 

To examine in more detail the coarsening we calcu- 
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FIG. 3: The time evolution of a random initial configuration 
of N = 5 ■ 10 s particles on a ladder of length L = 25 ■ 10 3 
and for g = 0.1. The configuration of every fifth particle is 
recorded every 250 MC steps. 



lated the average fraction of active particles ua = N^/N 
as a function of time for several values of q (Fig. 0|. 
Since active particles are mainly at the border of clus- 
ters [lH , their number Na is proportional to the number 
of clusters and the characteristic cluster size should scale 
as N/Na = On the semi-log plot of our data 

one can clearly see that 1/ua and thus the characteris- 
tic cluster size, after an initial fast growth, at late stage 
increases logarithmically slowly in time (~ log(t)). We 
also measured the average size of the maximum cluster 
N m ax just after the initial fast growth. Our data (in- 
set in Fig. 2]) suggest that this quantity diverges around 
q = q c ~ 0.041. 

Logarithmically slow coarsening suggests that the dy- 
namics of the model might exhibit some glassy features. 
Considering small- and large-g phases as analogs of high- 
and low-temperature phases, respectively, one can ex- 
amine evolution of our model under continuous cooling. 
In our simulations we increase q linearly in time from 
till 0.25 as q(t) — ct, where c can be interpreted as 
a cooling rate. The average value of ha as a function 
of q is shown in Fig. One can notice that for slow 
cooling riA rapidly decreases around q = 0.04 and such 
a value agrees with the divergence of maximum cluster 
size N max (inset in Fig. . In glassy dynamics it is of 
interest to examine the zero-temperature characteristic 
length lo and its dependence on the cooling rate. Some 
arguments [l4[ and Monte Carlo simulations [2(| sug- 
gest that in glassy systems lo should only very weakly 
increase with the decreasing cooling rate. Very small de- 
pendence on the cooling rate is also reported in some 
Molecular Dynamics simulations [2l[ as well as in some 
experiments [22j j . Linear fit to the average value of ua, as 
measured at the end of the cooling, as a function of the 
cooling rate c suggests that Ua(q — 0.25) ~ c 5 (inset 



FIG. 4: The time evolution of the inverse of the fraction of 
active particles ua (semi- log scale). Simulations were made 
for N = 2 ■ 10 4 and averaged over 10 2 independent runs. The 
inset shows the q dependence of the inverse of the maximum 
cluster size N ma x- These data suggest that N max diverges 
around q — q c ~ 0.041. 
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FIG. 5: The fraction of active particles ua as a function of 
q during the linear cooling with the values of cooling rate c 
indicated. Results are averages over 10 2 -10 4 independent runs 
(TV = 2 ■ 10 5 ). The inset shows the value of ua at the end of 
the cooling (q — 0.25) as a function of cooling rate c (log-log 
scale). The dotted line has a slope 0.5 that corresponds to 
n A ~ c - 5 . 



in Fig. [5]). Such a relation would imply that the charac- 
teristic cluster size at the end of the cooling increases as 

c -0.5_ 

Taking into account that the model has 
logarithmically-slow coarsening, such a fast increase of 
the characteristic length seems to be rather surprising. 
In our opinion, there are two possible scenarios that 



4 



might explain such a behaviour. First, there might 
be a regime close to the transition point q c < q < q' c 
where the coarsening dynamics is much faster (with 
presumably power-law coarsening). The dominant part 
of the coarsening would take place during the time 
spent in that regime and power-law increase of the 
characteristic length would be the expected feature. As 
shown by Shore et al., similar scenario takes place in a 
certain three-dimensional Ising model with next-nearest 
neighbour interactions [bij . Our results on coarsening in 
Fig. 2] extends only up to q = 0.07 and we cannot exclude 
that for smaller q a power-law coarsening would operate 
(i.e., q' c , if exists, should be smaller than 0.07). Another 
possibility is that the regime with logarithmically slow 
coarsening extends up to the transition point q = q c 
but the dynamics at this point being sufficiently fast to 
generate such a growth. Such a situation most likely 
occurs in three-dimensional gonihcdric Ising models and 
is caused by vanishing of energy barriers precisely at the 
critical temperature [13| . 

Although the dynamics of our model has much larger 
cooling-rate effects than that expected in glasses, we can- 
not exclude that some other one-dimensional system will 
be more realistic with this respect. Let us notice that 
in three-dimensional Ising model with 4-spin (plaquette) 
interaction, which is a particular version of the already 
mentioned gonihedric model, there are very small, and 
presumably logarithmic, cooling-rate effects (as well as 



logarithmically slow coarsening) [121 ]. 

In conclusions, using extensive numerical simulations 
we have shown that in the two-lane two-species driven 
diffusive particle model formation of clusters proceeds via 
logarithmically slow coarsening. Adopting phenomenol- 
ogy developed for glassy systems we implemented a 
cooling protocol for our model. Obtained results show 
that the model misses an important ingredient of glassy 
dynamics and despite slow coarsening it has rather 
strong, power-law cooling-rate effects. Together with 
some Molecular Dynamics simulations of off-lattice sys- 
tems [2l| , our work suggests that absence of glassy transi- 
tion might be a more general feature of one-dimensional 
systems. Of course, one cannot exclude that a certain 
modification of our model with e.g., different density of 
particles, different concentrations of "positive" and "neg- 
ative" particles, or jumping rules without CP symmetry 
will behave differently, exhibiting both slow coarsening 
and small cooling-rate effects. Alternatively, there might 
be some more fundamental reasons prohibiting glassy dy- 
namics in one-dimensional systems, and further research 
that would resolve that problem, in our opinion, should 
be undertaken. 
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